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Abstract 

Mathematical  modeling  has  been  extensively  applied  to  the  study  and  development  of  fuel  cells.  In  this  work,  the  objective  is  to  characterize  a 
mechanistic  model  for  the  anode  of  a  direct  ethanol  fuel  cell  and  perform  appropriate  simulations. 

The  software  Comsol  Multiphysics®  (and  the  Chemical  Engineering  Module)  was  used  in  this  work.  The  software  Comsol  Multiphysics®  is 
an  interactive  environment  for  modeling  scientific  and  engineering  applications  using  partial  differential  equations  (PDEs).  Based  on  the  finite 
element  method,  it  provides  speed  and  accuracy  for  several  applications. 

The  mechanistic  model  developed  here  can  supply  details  of  the  physical  system,  such  as  the  concentration  profiles  of  the  components  within 
the  anode  and  the  coverage  of  the  adsorbed  species  on  the  electrode  surface.  Also,  the  anode  overpotential-current  relationship  can  be  obtained. 
To  validate  the  anode  model  presented  in  this  paper,  experimental  data  obtained  with  a  single  fuel  cell  operating  with  an  ethanol  solution  at  the 
anode  were  used. 

©  2008  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  advantages  of  fuel  cells  as  clean,  silent  and  efficient 
power  sources  have  been  demonstrated,  and  proton  exchange 
membrane  fuel  cells  (PEMFCs)  are  presently  among  the  most 
promising  of  these  devices  [1,2].  Although  a  PEMFC  working 
with  hydrogen  presents  the  highest  efficiency  and  environmental 
compatibility,  hydrogen  is  not  the  best  choice  as  a  fuel  because 
of  the  problems  of  production,  storage  and  distribution.  A  liq¬ 
uid  fuel  like  a  low  molecular  weight  alcohol  can  be  used  in 
a  direct  alcohol  fuel  cell,  which  has  the  advantage  of  being 
a  much  simpler  system.  Methanol  has  been  used  as  fuel  with 
some  success,  but  it  has  some  disadvantages,  e.g.  it  is  relatively 
toxic,  inflammable  with  a  low  boiling  point  (65  °C),  and  it  is 
neither  a  primary  fuel,  nor  a  renewable  fuel.  Therefore,  other 
alcohols  are  being  considered  as  alternative  fuels.  Ethanol  is 
an  attractive  liquid  fuel  for  direct  alcohol  fuelled  systems.  It 
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is  the  major  renewable  biofuel  obtained  from  the  fermentation 
of  biomass,  and  ethanol  is  less  toxic  than  methanol.  Wang  et 
al.  [3]  compared  the  performance  of  fuel  cells  employing  an 
H3PO4 -doped  polybenzimidazole  membrane  operating  on  var¬ 
ious  methanol-alternative  fuels.  They  found  that  ethanol  is  a 
promising  alternative  fuel  with  an  electrochemical  activity  com¬ 
parable  to  that  of  methanol. 

Several  experimental  studies  on  the  electro-oxidation  of 
ethanol  have  been  devoted  mainly  to  identifying  the  adsorbed 
intermediates  on  the  electrode  and  elucidating  the  reaction 
mechanism  by  means  of  various  techniques,  as  differential 
electrochemical  mass  spectrometry  (DEMS),  in  situ  Fourier 
transform  infrared  spectroscopy  (FTIRS)  and  electrochemical 
thermal  desorption  mass  spectroscopy  (ECTDMS)  [4-9].  Based 
on  the  foregoing  work,  the  global  oxidation  mechanism  of 
ethanol  in  acid  solution  may  be  summarized  by  the  following 
scheme  of  parallel  reactions: 

CH3CH2OH  ->  [CH3CH2OH]ad  ^  Clad,  C2ad 

— >  C02  (total  oxidation)  (1) 
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Nomenclature 

Ck  species  concentration  (mol  m-3) 

£>k  diffusivity  (m2  s  - 1 ) 

F  Darcy’s  drag 

j  current  density  (A  m-2  or  mA  cm-2) 

jo  constant  in  Tafel  equation  (A  m-2) 

kO  kinetic  constant  (ms-1) 

kO'  kinetic  constant 

kO"  kinetic  constant  (mol  m-2  s-1) 

kl  kinetic  constant  (ms-1) 

kY  kinetic  constant 

M  kinetic  constant  (mol  m-2  s-1) 

k4'  kinetic  constant  (mol  m-2  s-1) 

k5  kinetic  constant  (mol  m-2  s-1) 

L  length  (m) 

R  resistance  (£2  m2)  =  0.285  x  10-4  Q  m2 
t  time  (s) 

u  component  of  the  velocity  vector  u  in  the  v  direc¬ 

tion  (ms-1) 

u  velocity  vector  (ms-1) 

v  component  of  the  velocity  vector  u  in  the  y  direc¬ 

tion  (ms-1) 
x  v  axis  (m) 

y  y  axis  (m) 

Greek  letters 

a  transfer  coefficient 

r  site  density  (mol  m-2) 

s  porosity 

rj  overpotential  (V  or  mV) 

0  coverages  of  the  different  species 

/x  viscosity  (kg  m- 1  s  - 1 ) 

p  density  (kg  m- 3 ) 


CH3CH2OH  ->  [CH3CH2OH]ad  ->  CH3CHO 

CH3COOH  (partial  oxidation)  (2) 

The  formation  of  C02  is  believed  to  go  through  two  adsorbed 
intermediates  Clad  and  C2ad,  which  represent  fragments  with 
one  and  two  carbon  atoms,  respectively.  In  spite  of  many 
advances  in  the  understanding  of  the  mechanism  of  ethanol  oxi¬ 
dation,  there  are  still  some  unclear  aspects.  For  instance,  there  is 
some  controversy  on  whether  acetic  acid  is  formed  in  one  step  or 
through  the  aldehyde.  Also,  there  is  no  agreement  regarding  the 
nature  of  the  adsorbed  species.  According  to  some  workers,  the 
carbon-carbon  bond  is  preserved,  so  a  larger  quantity  of  inter¬ 
mediates  of  the  type  C2  are  formed  [4,5],  but  others  claim  that 
the  main  intermediates  contain  only  one  carbon  atom  and  are  of 
the  type  Cl  [6,8].  Breaking  the  C-C  bond  for  a  total  oxidation  to 
C02  is  a  major  problem  in  ethanol  electrocatalysis.  Thus,  high 
yields  of  partial  oxidation  products,  CH3CHO  and  CH3COOH, 
are  found  at  Pt  catalysts  [10,11].  These  parallel  reactions  cause 
a  considerable  lowering  of  the  fuel  capacity  to  generate  electric¬ 


ity  and  produce  undesirable  substances.  Platinum  is  at  present 
the  best-known  catalyst  for  the  adsorption  and  dissociation  of 
small  organic  molecules,  but  it  has  a  limited  ability  for  break¬ 
ing  the  C-C  bond.  In  this  context,  the  research  objectives  in 
this  area  have  been  the  search  for  more  effective  electrocata¬ 
lysts,  which  ideally  promote  the  complete  oxidation  of  ethanol 
at  potentials  as  low  as  possible.  Presently,  research  in  electro¬ 
catalysts  is  focused  on  the  development  of  binary  and  ternary 
Pt-based  materials. 

Mathematical  modeling  has  been  extensively  applied  in  the 
study  and  development  of  fuel  cells  [12].  Modeling  studies  of 
different  physical  systems,  e.g.  usual  cathodes  and  biochemical 
fuel  cells  [13],  have  been  developed  in  this  laboratory.  In  this 
work,  a  mechanistic  model  for  an  ethanol  anode  is  proposed  and 
appropriate  simulations  performed,  with  the  objective  of  provid¬ 
ing  simulation  elements  that  can  help  with  the  understanding  of 
the  anodic  process  in  direct  ethanol  fuel  cells. 

2.  Methodology 

Simulations  were  performed  using  the  software  Comsol 
Multiphysics®  [14].  This  software  is  an  interactive  environment 
for  modeling  scientific  and  engineering  systems  based  on  partial 
differential  equations  (PDEs).  Using  the  finite  element  method 
[15,16],  it  provides  speed  and  accuracy  for  multiphysics  appli¬ 
cations  and  it  covers  all  the  facets  of  the  modeling  process.  It 
contains  CAD  (“Computer  Aided  Design”)  tools,  interfaces  for 
specifications  of  equations  and  physical  and  chemical  parame¬ 
ters,  mesh  generation,  optimized  solvers,  as  well  as  visualization 
and  post-processing  tools. 

2.7.  Specialized  modules 

This  work  was  developed  using  the  Chemical  Engineering 
Module  [17],  which  is  divided  into  many  application  modes. 

2.7.7.  Momentum  balances 

The  Navier-Stokes  equations  for  fluid  flow; 

Darcy’s  Law  for  porous  media  flow. 

2.7.2.  Energy  balances 

Heat  transfer  by  diffusion,  convection  and  conduction. 

2.1.3.  Material  balances 
Convection  and  diffusion. 

In  these  applications,  the  physical  (and  chemical)  proper¬ 
ties  appearing  in  the  equations  are  specified  so  the  PDEs  are 
formulated. 

The  modeling  procedure  using  the  Chemical  Engineering 
Module  is  analogous  to  the  general  methodology  of  COMSOL 
Multiphysics®  [18].  It  involves  the  following  steps: 

1.  create  or  import  the  geometry  in  ID,  2D  or  3D; 

2.  select  the  equations  that  define  the  system; 

3.  specify  the  properties  or  coefficients  in  the  selected  equa¬ 
tions; 

4.  set  the  boundary  and  initial  conditions; 
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5.  generate  and  refine  the  finite  element  mesh; 

6.  run  the  simulation; 

7.  visualize  and  process  the  results. 

For  the  validation  experiments,  the  membrane  electrode 
assembly  (MEA)  was  prepared  by  hot  pressing  a  pretreated 
Nation®  117  membrane  between  a  Pt/C  (metal  loading  30  wt.%) 
E-TEK  cathode  (1  mg  cm-2  Pt)  and  a  Pto.ssSno.sWo.v/C  (metal 
loading  40  wt.%)  anode  (lmgcm_2Pt),  at  125  °C  for  2  min 
under  a  pressure  of  35  kg  cm-2.  The  operating  fuel  cell  perfor¬ 
mance  was  determined  in  a  single  direct  ethanol  fuel  cell  (DEFC) 
with  an  electrode  geometric  surface  area  of  4.6  cm2.  The  tem¬ 
perature  was  set  to  90  °C  for  the  fuel  cell  and  95  °C  for  the 
oxygen  humidifier.  The  distribution  of  reaction  products  under 
the  DEFC  operating  conditions  was  determined  at  a  current  of 
30  mA  cm-2.  The  final  products  were  trapped  on  the  anode  and 
the  products  derived  from  the  crossover  process  were  collected 
on  the  cathode  outlet  in  an  iced  flask.  A  High  Performance  Liq¬ 
uid  Chromatograph  (HPLC)  using  an  Aminex  HPX-87H  column 
from  Bio-Rad  with  a  UV  detector  SPD-MICA  from  Shimadzu 
was  used  to  analyze  quantitatively  the  reaction  products  at  the 
outlet  of  both  anode  and  cathode. 

3.  Results  and  discussion 

Fig.  1  shows  the  simplified  scheme  of  the  fuel  cell  anode  used 
here.  This  2D  pattern  is  usually  used  to  simulate  parallel  flow 
fields  (interdigitated  flow  fields,  on  the  other  hand,  are  usually 
associated  with  3D  models). 

The  mathematical  model  of  the  ethanol  anode  is  based  on 
material  balances  considering  diffusive  and/or  convective  fluxes 
((3(Ck))/30  =  —  V(wCk)  +  V(eDkVCk)  =  0  for  the  species  in 
solution  (mainly  water,  ethanol,  acetic  acid  and  acetaldehyde 
or  CO2).  Additionally,  the  model  assumes  Fickis  Law  for  dif¬ 
fusion  and  Tafel  kinetics  coupled  to  adsorption  steps  on  the 
catalytic  surface.  The  treatment  is  similar  to  that  of  Nordlund 
and  Lindberg  [19]  for  a  methanol  anode.  Besides  Nordlund  and 
Lindberg  [19],  it  is  relevant  to  mention  the  works  of  Zhou  et 
al.  [20],  Sundmacher  et  al.  [21],  and  Vidakovic  et  al.  [22].  The 
velocity  field  was  also  calculated  (according  to  ((3(pw))/3f)  = 
—V(puu)  —  VP  +  V(/x(Vm  +  (Vw)T))  +  F). 

Reaction  steps  on  the  catalytic  surface  are  the  scheme  of  Eqs. 
(3a)-(3e).  Eq.  (3d)  represents  the  formation  of  the  species  OHads, 


y 


catalyst  layer 


porous  layer 


channel 


x 


which  plays  a  very  important  role  in  the  oxidation  of  some  inter¬ 
mediate  species  by  providing  the  necessary  oxygen  atoms.  On 
the  real  catalyst  surface  of  binary  and  ternary  Pt-based  electro¬ 
catalysts,  there  are  different  sites  with  different  crystal  structures 
available  for  adsorption  and  reaction.  Small  organic  molecules 
will  mostly  adsorb  on  Pt.  Different  sites  may  present  different 
rate  constants,  but  very  detailed  information,  not  available  at  the 
moment,  is  necessary  to  differentiate  them  quantitatively.  So, 
we  have  chosen  as  a  good  initial  approach  to  describe  the  cata¬ 
lyst  through  one  set  of  rate  constants,  such  as  in  Nordlund  and 


Lindberg  [19]. 

CH3CH2OH  CH3CHOads  +  2H+  +  2e“  k().k(Y  (3a) 

CH3CHOads^  CH3CHO  /(f)"  (3b) 

CH3CH2OH  +  H20 

CH3COads  +  2H+  +  3e“+H30+  kl,kl'  (3c) 

H20  OHads  +  H++e“  M,M'  (3d) 

CH3COads  +  OHads  -*  CH3COOH  k5  (3e) 


For  the  coverages  0  of  the  different  species  on  the  catalyst  sur¬ 
face,  the  following  steady-state  conditions  are  postulated: 

r—  dfH°  =  k0 '  Cch''ch^oh  ‘ (1  _  m  ‘  exP(0-5  ’ const  r]} 

-  k0'  •  6»ch3cho  •  exp(— 0.5  •  const,  rj) 

-  kO"  ■  0ch3cho  =  0  (4a) 

rCWc^CO  =  kl  •  CCH3ch2oh  •  (1  -  Z0)  •  exp(al  •  const,  r,) 

-  k Y  ■  0ch3co  ■  exp(— (1  -  al)  •  const.  rj) 

-  k 5  ■  0oh  ■  #ch3co  =  0  (4b) 

d$OH 

L  — —  =  M  •  (1  -  Z6)  •  exp(0.5  •  const,  rj) 
d  t 

—  k4f  •  $oh  •  exp(— 0.5  •  const,  r]) 

—  k5  •  0ou  •  $ch3co  =  0  (4c) 

From  which  the  variables  0  are  obtained  as  a  function  of  the 
anode  overpotential. 

The  scheme  corresponding  to  Eqs.  (3a)-(3e)  is  adequate  for 
the  electro-oxidation  of  ethanol  on  a  platinum  electrode,  where 
the  main  pathways  lead  to  the  formation  of  acetaldehyde  and 
acetic  acid,  so  the  scheme  corresponds  to  a  realistic  model. 

Because  modeling  allows  to  postulate  useful,  albeit  unre¬ 
alistic  pathways  for  the  reaction,  an  ideal,  unrealistic  model 
consisting  in  the  total  oxidation  process  will  also  be  considered 
as  reference: 

CH3CH2OH  +  3H20  2C02  +  12H+  +  12e“  (5) 

As  boundaries,  the  following  conditions  were  considered:  for 
0<y< channel  height:  velocity  components  and  species  con- 


Fig.  1.  Scheme  of  a  porous  anode. 
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centration  at  channel  entrance: 


for  the  catalyst  interface: 


U\x= 0  =  Win,  V\x=0  =  0 
Ch2oIx=o  =  50111  mol  m-3, 

Cch3ch2ohIx=o  =  2000  mol  m-3, 

Cqo2  U=o  =  0  (no  CO2  at  channel  entrance) 

or 

Ch2oL=o  =  50111  mol  m-3, 

Cch3ch2ohL=o  =  2000  mol  m-3, 

^acetaldehyde I x=o  =  0  (no  acetaldehyde  at  channel  entrance), 
^acetic  atidL=o  =  0  (no  acetic  acid  at  channel  entrance) 

for  channel  height  <  y  <  total  height: 


-sDk  («*»)  I y= total  height  =  diffusive  flux  = 

=  reaction  fluxspecie, 

(  9(Cch3ch2oh)\  .  .  a 

sDk  I  - — - j  v=total  height  =  diffusive  flux 

=  JjJ  =  reaction  fluxspeciei 

~sDk  )  j  l  v=totai  height  =  diffusive  flux  =  -2- 

=  reaction  fluxspecie, 

Cch3ch2oh  \0,25  ( a  A  Fy\ 

CcH3CH2OH,ref/  \  RT  ) 

diffusive  fluxspecie.  =  adsorption/reaction  fluxspecie- 


u\x=o  =  0,  v\x=o  —  0  (zero  velocity  at  “wall”), 
^  |x_0  =  0  (insulation/ symmetry) 
u\x=l  —  0,  v\ x=l  =  0  (zero  velocity  at  “wall”), 
^  \X=L  =  0  (insulation/ symmetry) 

for  0  <  v  <  cell  length: 


W  | y=0  =  0, 


v\y= o  =  0  (zero  velocity  at  “wall”), 
=  0  (insulation/ symmetry) 


The  operating  conditions  and  the  parameters  considered  are 
presented  in  Table  1 . 

A  finite  element  mesh  of  75  x  17  elements  was  used  in  the 
simulations  when  considering  the  ideal  model  and  a  mesh  of 
80  x  17  elements  was  used  when  considering  the  realistic  model. 
In  the  v  “along-channel”  direction  (see  Fig.  1),  around  100  grid- 
points  are  typically  required.  We  used  75  (ideal  model)  and  80 
(realistic  model).  Besides,  6-10  grid-points  are  needed  to  suffi¬ 
ciently  resolve  each  of  the  distinctive  regions  “through-plane”  (y 
direction;  channel,  diffusion,  membrane. . .).  We  used  10  (chan¬ 
nel)  +  7  (diffusion).  So,  we  used  very  typical  values.  Considering 
the  reference  model,  that  is  the  total  oxidation  process,  Figs.  2-4 
show  the  concentration  profiles  for  the  species  in  solution  at 
an  anode  overpotential  of  0.37  V  and  the  conditions  given  in 


Table  1 

Parameters  used  in  the  modeling 


Water  diffusivity  (in  liquid  phase) 

T>h2o 

5  x  10-10  m2  s-1 

Ref.  [23] 

Ethanol  diffusivity  (in  liquid  phase) 

A:h3ch2oh 

6  x  10-9  m2  s-1 

Ref.  [23] 

C02  diffusivity  (in  liquid  phase) 

Dco2 

8  x  10-9  m2  s-1 

Ref.  [23] 

Acetaldehyde  diffusivity  (in  liquid  phase) 

£>ch3cho 

6  x  10-9  m2  s-1 

Assumed  (equal  to  DCh3cooh) 

Acetic  acid  diffusivity  (in  liquid  phase) 

A:h3cooh 

6  x  10-9  m2  s-1 

Ref.  [23] 

Diffusion  layer  porosity 

£ 

0.9  (typical  value)  or  1  (in  the  simulations 
considering  the  ideal  model) 

Tafel  equation  constant  (ethanol,  ideal  model) 

jo 

68.33  x  10-1  Am“2 

Ref.  [24] 

Ethanol  reference  concentration  (ideal  model) 

CcH3CH2OH,ref 

500  mol  m-3 

Ref.  [24] 

Anode  transfer  coefficient  (ethanol,  ideal  model) 

Ol  A 

0.48 

Ref.  [24] 

Kinetic  constant 

kO 

3  x  10-6  ms-1 

Assumed 

Kinetic  constant 

k0' 

0 

Reaction  (3a)  considered  irreversible 

Kinetic  constant 

kO" 

1  x  10-2  molm-2  s_1 

Ref.  [25] 

Kinetic  constant 

kl 

7  x  10-8  ms-1 

Assumed 

Kinetic  constant 

kV 

0 

Reaction  (3c)  considered  irreversible 

Coefficient 

al 

0.8 

Assumed 

Kinetic  constant 

k\ 

1  x  10-5  molm-2  s_1 

Ref.  [22] 

Kinetic  constant 

M' 

1  mol  m-2  s_1 

Ref.  [22] 

Kinetic  constant 

k5 

1  x  10_1  molm-2  s_1 

Assumed 

Channel  height 

0.0008  m 

Measured 

Diffusion  layer  height 

0.0004  m 

Typical  value 

Anode  length 

L 

0.0215  m 

Measured 

Temperature 

T 

90  °C  (363  K) 

Measured 
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Contour:  Concentration,  cl  [md/m3] 


x  (m) 


Max:  4.982e4 

xlO4 

^4.982 

- 4,925 

- 4,868 

- 4.81 

- 4,753 

-  A  696 

- 4.638 

- 4.581 

- 4.524 

-  4,466 

4.409 

4,352 

4.294 

4.237 

4.18 

4.122 

4.065 

4.008 

3,95 

- 3.893 

Min:  3.893e4 


Fig.  2.  Water  concentration  profile  for  the  reference  model. 


Table  1.  Water  consumption  can  be  seen  in  Fig.  2,  but  the  con¬ 
centration  does  not  become  too  low  in  any  place.  The  shape  of 
the  curves  reveals  the  effect  of  a  predominantly  convective  flux 
in  the  channels  and  a  mainly  diffusive  flux  in  the  porous  layer. 
Finally,  Fig.  5  shows  the  overpotential-current  curve  that  results 
for  the  ideal  model. 

Considering  now  the  realistic  model,  Figs.  6-9  show  the  con¬ 
centration  profiles  at  an  anode  overpotential  of  0.37  V.  When 
comparing  Figs.  2  and  6,  it  is  observed  that  water  concentration 
close  to  catalyst  surface  is  larger  for  the  realistic  model.  This  is 
correct,  because  in  the  real  situation  less  water  is  consumed.  The 
ideal  model  predicts  higher  current  densities,  which  corresponds 
to  higher  water  consumption.  When  comparing  Figs.  3  and  7,  it 
is  observed  almost  the  same  consumption  of  ethanol,  but  as  will 


be  seen  in  Fig.  11,  the  current  density  for  the  realistic  model  will 
be  lower. 

Fig.  10  shows  the  coverages  of  the  adsorbed  intermedi¬ 
ate  species  postulated  in  Eqs.  (3a)-(3e),  as  a  function  of  the 
overpotential.  The  results  show  that  at  an  overpotential  around 
0.41V,  the  surface  coverage  of  OH  starts  being  larger  than 
the  surface  coverage  of  CH3CO.  Finally,  Fig.  11  shows  the 
overpotential-current  plot  for  the  realistic  model  compared 
with  that  obtained  for  the  ideal  model.  It  can  be  seen  that  at 
higher  overpotentials  (over  0.44  V),  the  increasing  OH  coverage 
decreases  the  cell  performance.  When  this  non-Tafel  behaviour 
is  observed,  it  means  that  adsorbed  hydroxide  can  compete  with 
the  adsorption  of  other  species.  Nevertheless,  it  is  important 
to  notice  that  a  model  which  describes  the  catalytic  reaction 


Fig.  3.  Ethanol  concentration  profile  for  the  reference  model. 


288 


R.  Sousa  Jr  et  al.  /  Journal  of  Power  Sources  180  (2008)  283-293 


Contour:  Concentration,  c3  [md/m3] 


x(m) 
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Fig.  4.  CO2  concentration  profile  for  the  reference  model. 


through  one  set  of  rate  constants  (adsorption  on  the  same  type 
of  site),  such  as  in  this  paper  and  in  Nordlund  and  Lindberg 
[19],  that  kind  of  behaviour  (competition)  can  be  “favoured”. 
So,  if  the  non-Tafel  behaviour  is  not  observed,  a  more  rig¬ 
orous  model,  with  a  more  complex  description  of  adsorption 
sites  can  improve  the  modeling  predictions  at  high  overpoten¬ 
tials.  However,  as  mentioned  above,  very  detailed  information 
is  necessary  to  differentiate  rate  constants  for  different 
sites. 

To  validate  the  anode  model  presented  in  this  paper,  exper¬ 
imental  data  from  this  laboratory  obtained  with  a  single  cell 
operating  with  an  ethanol  anode  and  a  serpentine  flow  field 
were  used.  Since  a  serpentine  flow  field  was  used,  we  consid¬ 
ered  the  scheme  in  Fig.  1  to  model  the  flow  path  as  a  sequence 
of  four  flow  steps  (three  U-turn  regions  not  affecting  the  con¬ 
centration  profiles,  only  changing  the  flow  direction),  extending 


Fig.  5.  Overpotential-current  plot  for  the  reference  model. 


the  modeling  approach.  An  oxygen  cathode  (at  3  atm)  and  a 
Nation®  117  membrane  completed  the  system  (see  Fig.  12). 
It  is  important  to  notice  that  the  model  is  still  2D,  with  one 
v  “along-channel”  dimension,  one  y  “through-plane”  dimen¬ 
sion  and  no  dimension  in  the  direction  perpendicular  to  the 
xy  plane,  because  the  effect  of  the  U-turn  on  the  flow  and 
concentration  pattern  was  not  considered  [26].  That  seems  to 
be  the  only  possible  way  of  considering  a  2D  approach  for 
the  problem,  because  the  along-channel  x  direction  and  the 
through-plane  y  direction  are  the  more  important  and  can¬ 
not  be  omitted.  However,  a  complete  modeling  approach  must 
always  be  3D  (but  3D  models  require  additional  computer 
capabilities). 

A  standard  mathematical  model  of  a  cathode  was  considered. 
It  is  based  on  species  balance  equations.  Other  considerations 
are:  steady-state  operation,  Fick’s  diffusion  mechanism  and 
Tafel  kinetics.  Regarding  the  Nation®  membrane,  it  was  mod¬ 
eled  as  permeable,  in  a  sense  that  crossover  is  permitted.  Of 
course,  the  inclusion  of  the  membrane  and  the  crossover  pro¬ 
cess  will  result  in  a  distribution  of  species  not  only  within  the 
anode,  but  also  in  the  whole  cell  and  particularly  at  the  cathode. 
The  standard  parameters  considered  to  model  the  cathode  are 
presented  in  Table  2. 

The  experiments  with  the  single  cell  were  done  at  a  current 
density  of  300  A  m-2  (30  mA  cm-2). 


Table  2 

Parameters  considered  to  model  the  cathode 


Oxygen  diffusivity 

Dq2 

0.08  x  10"4  m2 

Diffusion  layer  porosity 

£ 

0.9 

Constant  in  Tafel  equation 

jo/ Co2,  ref 

0.0064  A  mmol 

Cathode  transfer  coefficient 

ac 

2 

Channel  height 

0.0008  m 

Gas  diffusion  layer  height 

0.0004  m 

Cathode  length 

L 

0.0215  m 
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Max:  4.981e4 

xlO4 
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Fig.  6.  Water  concentration  profile  for  the  realistic  model. 


The  results  in  terms  of  the  concentrations  of  the 
product  species  (ethanol  concentration  was  not  measured) 
were 

at  the  ANODE  channel  exit: 
acetic  acid  =  7  mol  m-3 
acetaldehyde  =  8  mol  m-3, 

at  the  CATHODE  channel  exit: 
acetic  acid  =  59  mol  m-3 
acetaldehyde  =  31  mol  m-3 . 

Figs.  13a,  b  and  14a,  b  show  the  simulation  results  at  a  current 
density  of  290  Am-2  (29  mA  cm-2),  in  terms  of  the  concentra¬ 


tion  profiles  of  acetic  acid  and  acetaldehyde  at  the  anode  and 
cathode  exits. 

The  simulated  average  acetic  acid  and  acetaldehyde  con¬ 
centrations  at  the  anode  channel  exit  were  6  and  9molm-3, 
respectively.  Average  acetic  acid  and  acetaldehyde  concentra¬ 
tions  at  the  cathode  channel  exit  were  20  and  66  mol  m-3, 
respectively.  So,  it  may  be  concluded  that  there  is  a 
good  agreement  between  the  experimental  and  the  sim¬ 
ulated  results  for  the  current  density,  the  cell  potential 
(V=  0.85  —  ?7anode  +  ^cathode  ~  Rj)  and  the  concentrations  of 
product  species  at  the  anode.  On  the  other  hand,  there  is 
obviously  some  disagreement  regarding  the  concentrations  of 
product  species  at  the  cathode.  The  experimental  results  show 
that  acetaldehyde  concentration  at  the  cathode  is  smaller  than 


Contour:  Concentration,  c2  [mol/m3] 
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Fig.  7.  Ethanol  concentration  profile  for  the  realistic  model. 
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Fig.  8.  Acetic  acid  concentration  profile  for  the  realistic  model. 


Contour:  Concentration,  c5  [mol/m3] 
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Fig.  9.  Acetaldehyde  concentration  profile  for  the  realistic  model. 


Fig.  10.  Coverages  of  the  adsorbed  species  on  the  catalyst  surface  as  a  function  Fig.  11.  Overpotential-current  plot  for  the  realistic  model  compared  with  the 
of  the  anode  overpotential  for  the  realistic  model.  reference  model. 
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Fig.  12.  Scheme  of  a  unitary  cell. 


the  simulated  value,  while  the  concentration  of  acetic  acid  is 
larger  than  the  simulated  value.  This  can  be  understood  if  it 
is  considered  that  in  the  experiments  some  acetaldehyde  may 
have  been  lost  by  evaporation,  because  of  its  high  volatility,  and 
also  that  some  acetaldehyde  may  have  been  oxidized  to  acetic 
acid. 

Finally,  experimental  data  from  this  laboratory  (fuel  cell 
potential  vs.  current  density,  with  the  same  ternary  anode)  were 
used  to  carry  out  a  comparison  between  experimental  and  model¬ 
ing  results  for  a  complete  polarization  curve,  as  shown  in  Fig.  15. 
Again,  a  good  agreement  between  model  predictions  and  exper¬ 
imental  results  was  obtained.  It  is  important  to  mention  that  the 
comparison  of  modeling  with  experimental  results  for  the  fuel 
cell  is  in  the  operational  (anode)  overpotential  range,  which  is 
below  0.44  V  (and  below  the  non-Tafel  behaviour  region  men¬ 
tioned  above).  Fig.  16  shows  separately  the  anode  and  cathode 
overpotentials. 


Fig.  13.  (a)  Simulation  of  a  unitary  cell  operating  with  an  ethanol  anode  (acetic  acid  concentration  profile — cathode  exit),  (b)  Simulation  of  a  unitary  cell  operating 
with  an  ethanol  anode  (acetic  acid  concentration  profile — anode  exit). 
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Fig.  14.  (a)  Simulation  of  a  unitary  cell  operating  with  an  ethanol  anode  (acetaldehyde  concentration  profile — cathode  exit),  (b)  Simulation  of  a  unitary  cell  operating 
with  an  ethanol  anode  (acetaldehyde  concentration  profile — anode  exit). 
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Fig.  15.  Fuel  cell  potential  vs.  current  density. 


Fig.  16.  Anode  and  cathode  overpotentials. 
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4.  Conclusions 

A  mechanistic  model  for  the  anode  of  a  direct  ethanol  fuel 
cell  was  developed  on  the  basis  of  a  realistic  mechanism  for  the 
reaction  that  considers  the  formation  of  acetaldehyde  and  acetic 
acid  among  the  products.  From  the  model,  it  was  possible  to  get 
details  of  the  physical  system,  such  as  the  concentration  profiles 
of  the  components  within  the  diffusion  anode  and  the  coverage 
of  the  adsorbed  species  on  the  electrode  surface. 

The  model  allowed  to  obtain  the  overpotential  vs.  current 
plot.  At  all  overpotentials,  the  currents  are  lower  than  those  that 
would  be  obtained  with  an  ideal  model  that  considers  the  total 
oxidation  to  CO2. 

The  concentrations  of  acetaldehyde  and  acetic  acid  at  the 
anode,  due  to  the  oxidation  reaction,  and  at  the  cathode,  due 
to  crossover,  showed  good  agreement  with  those  determined 
experimentally  in  a  single  direct  ethanol  fuel  cell.  Finally, 
the  comparison  between  experimental  data  from  this  labo¬ 
ratory  (fuel  cell  potential  vs.  current  density)  and  modeling 
results  showed  good  agreement,  giving  support  to  the  procedures 
employed  here. 
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